Parton rescattering and screening in Au+Au collisions at RHIC 
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We study the microscopic dynamics of quarks and gluons in relativistic heavy ion coUisions in the 
framework of the Parton Cascade Model. We use lowest order perturbative QCD cross sections with 
fixed lower momentum cutoff p™'". We calculate the time-evolution of the Debye-screening mass 
for Au+Au collisions at y/s — 200 GeV per nucleon pair. The screening mass is used to determine 
a lower limit for the allowed range of p™'". We also determine the energy density reached through 
hard and semi-hard processes at RHIC, obtain a lower bound for the rapidity density of charged 
hadrons produced by semihard interactions, and analyze the extent of perturbative rescattering 
among partons. 



The first experimental results from tlie Relativistic 
Heavy Ion Collider (RHIC) have generated a vast amount 
of interesting data |^ . A variety of theoretical models 
has been invoked to describe the observed phenomena, 
e.g., fluid dynamical models, perturbative QCD scatter- 
ing models, as well as models based on parton satura- 
tion and statistical approaches. Although these mod- 
els, which all contain adjustable parameters, have been 
fairly successful within their respective regimes of antici- 
pated applicability, they all have certain limitations. For 
example, fluid dynamics cannot describe transport phe- 
nomena occurring prior to local equilibration of the pro- 
duced matter, and it must fail above a certain, though 
unknown, value of px- Perturbative parton scattering 
models cannot describe the physics of equilibration and 
the formation of collective flow. Initial state parton sat- 
uration models do not include the final state interaction 
among partons, which leads to equilibration. 

The parton cascade model (PCM) was proposed, about 
a decade ago, with the aim to provide a unified frame- 
work for the description of phenomena involving high and 
intermediate values of pt • It is based on the premise 
that the microscopic dynamics during the early stages 
of a high-energy nuclear collision can be described as 
a cascade of two-body interactions among perturbative 
quarks and gluons, with certain modifications caused by 
the presence of a hot and dense medium. Because of 
its perturbative nature, the PCM cannot describe the 
hadronization stage, unless additional phenomenological 
models are introduced. 

Since the original formulation and implementation of 
the PCM by Geiger (VNI) Q, significant new insights 
into the dynamics of dense systems of partons have been 
gained. First and foremost among these is the applicabil- 
ity of semiclassical methods for the description of the par- 
tonic structure in the initial state and the earliest phase 
of its evolution ||^ . Another relevant development is the 
recognition of the importance of radiative processes for 
the attainment of local thermal equilibrium in the rapidly 



expanding matter. Radiative processes were already in- 
cluded in the original PCM in the leading-logarithmic 
approximation (LLA). The concept of initial-state sat- 
uration of the parton distribution, however, provides a 
novel feature which can be utilized to address the in- 
frared divergences of the perturbative PCM. 

In this letter we focus on the early, pre-thermal reac- 
tion phase of a heavy-ion reaction and address the fol- 
lowing questions: 

• What limits on the cut-off p^™, required by the 
infra-red divergence of the pQCD cross section, are 
imposed by the internal consistency of the model 
and by existing data on particle distributions? 

• What energy-density and multiplicity is reached 
through hard and semi-hard processes (alone) at 
RHIC within these limits? 

• Do partons undergo multiple perturbative binary 
collisions? Are these sufficient to produce (measur- 
able) collective effects? 

We will here restrict our discussion of the results from 
the Parton Cascade Model implementation (VNI/BMS, 
an improved and corrected version of the VNI implemen- 
tation [§)) to binary processes in leading order pQCD - a 
more detailed calculation involving NLO corrections such 
as radiative processes will be presented in a forthcoming 
publication 

The fundamental assumption underlying the PCM is 
that the state of the dense partonic system can be 
characterized by a set of one-body distribution func- 
tions Fi{x^,p"). Here i denotes the fiavor index (i — 
g,u,u,d,d, . . .) and x^^p" are coordinates in the eight- 
dimensional phase space. The partons are assumed to 
be on their mass shell, except before the first scatter- 
ing. In the numerical implementation, the continuous 
parton distribution functions (here we choose the GRV- 
HO parametrization |^]), are represented by (test) parti- 
cles 
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where (r) and gf (r), respectively, denote the space- 
time position and four-momentum of particle i. r is a 
variable (proper time) parameterizing the world-line of a 
particle. The factor = 0, 1 allows for the creation and 
annihilation of partons. 

Partons generally propagate on-shell and on straight- 
line paths between scattering events. Before their first 
collision, partons may have a space-like four-momentum, 
especially if they are assigned an "intrinsic" transverse 
momentum. 

The time-evolution of the parton distribution is gov- 
erned by a relativistic Boltzmann equation: 



with q and q' denoting different quark flavors. The am- 
plitudes for these processes have been calculated in refs. 
[H, ^ for massless quarks. The corresponding scattering 
cross sections are expressed in terms of spin- and colour- 
averaged amplitudes \M\'^: 
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For the transport calculation we also need the total cross 
section as a function of s which can be obtained from (j^): 



(2) 



where the collision term d is a nonlinear functional of 
the phase-space distribution function. Although the col- 
lision term, in principle, includes factors encoding the 
Bose-Einstein or Fermi-Dirac statistics of the partons, 
we neglect those effects here. 

The collision integrals have the form: 

(3) 



with 



and 



D{Fk{x,p)) - \{Fk{x,p) - \{Fk{x,p) (4) 
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Si is a statistical factor defined as S", = J] K''^^-Kl'^^\ 

with K™'°^'^ identical partons of species a in the initial 
or final state of the process, excluding the zth parton. 

The matrix elements \M.^ account for the following 
processes: 



(6) 



aa - 


aa 


aa - 


-> qq 


qa - 


qa 


qq' - 


Qq' 


qq - 


-> qq 


qq- 


- q'q' 


qq~ 


qq 


qq- 


aa 






qa - 


97 


qq~ 




qq- 





c.d 



da 



(8) 



ab—>-cd 



The low momentum-transfer cut-off p^'"" is needed to 
regularize the IR-divergence of the pQCD parton-parton 
cross section. A novel feature of our treatment involves 
the introduction of a factor /(x, Qq) for every 

primary parton involved in a binary collision in expres- 
sion ^ to account for the difference between the initializa- 
tion scale Qq and the scattering scale Q. A more detailed 
description of our implementation is in preparation |^ . 

One of the most crucial parameters of the PCM is the 
low momentum transfer cut-off p™™. Under certain as- 
sumptions this parameter can be determined from ex- 
perimental data for elementary hadron-hadron collisions 
]lO| , pr] , p^ . In the environment of a heavy-ion collision, 
colour screening will destroy the association of partons to 
particular hadrons, since for a sufficiently high density of 
colour charges, the colour screening radius becomes much 
smaller than the typical hadronic scale. It is therefore by 
no means clear whether the p™™ values extracted from 
hadron-hadron collisions are applicable to heavy-ion col- 
lisions. Calculating the colour screening mass fiD for a 
set of systems where interactions are governed by specific 
values of p™™ may allow us to determine a lower bound- 
ary for the allowed range of p™™ values, since only values 
of p™™ > fJ-D are physical. 

Following ref. , we use perturbative QCD to obtain 
the time evolution of the screening mass (t) ■ The par- 
ton cascade model provides the phase space distribution 
of the partons. The general form for the colour screening 
mass in the one loop approximation is |l^ |l5| 
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where as is the strong coupling constant, the Fi spec- 
ify the phase space density of gluons, quarks, and anti- 
quarks and q runs over the flavour of quarks. It is easy to 
verify that in the case of an ideal gas of massless partons. 



where the Fi reduce to Bose-Einstein or Fermi-Dirac dis- 
tributions (with vanishing baryochemical potential /is), 
Eq. ^ reduces to the standard result for the thermal De- 
bye mass [isl M. The partonic distribution will be ini- 
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FIG. 1: Top: time-evolution of the screening mass fin for 
several values of p™*" in Au+Au collisions at RHIC. Bottom: 
fj,D evaluated at r = l/p™'" as a function of p™'". The 
gray shaded area symbolizes the unphysical region in which 
prpin < g^jjjj^ therefore serves as a boundary for the allowed 
range of p™'" in Au+Au collisions at RHIC. 



tially anisotropic with respect to the beam axis and thus 
the screening mass of a gluon in the matter may depend 
on the direction of propagation. We have found that 

fxj^ differs from by 10% at most and we have there- 
fore assumed i^d — in the following discussion. We 
also note that the assumptions underlying this method 
are not strictly applicable to very early times, r < Az, 
where Az is the Lorentz contracted width of the nuclei. 



We begin by calculating the screening mass for par- 
tonic matter in central Au+Au collisions at RHIC {^/s = 
200 GeV per nucleon pair) . The upper frame of figure |l| 
shows the time evolution of the (transverse) screening 
mass, calculated according to eq. for several different 
values of p™'". The time t — Q marks the maximum 
overlap of the two colliding nuclei. The time-evolution of 
clearly reflects the dynamics of the collision: at early 
times the density of the system is large, leading to a large 
value the screening mass ^d- For later times, density and 
collision rate (see also figures || and ^ decrease, reducing 
the value of /iu by more than a factor of two. The strong 
time-dependence of the screening mass during the early 
pre-equilibrium phase will have an influence on the equi- 
libration process that is beyond the scope of our present 
work. A decreasing screening mass implies a rising cross 
section and thus enhanced multiple rescattering. This 
aspect is missing in many implementations of the PCM, 
which assume a time-independent screening mass for the 
entire duration of the collision |l^, |l^, ^ . 

A self-consistent calculation would utilize the time- 
dependent Debye-mass as the regulator of the interaction 
among partons instead of introducing a flxed cut-off p™*" 
as used in our present work. This would be achieved by 
calculating the parton-parton cross sections in the frame- 
work of the hard thermal loop (HTL) approximation , 
which accounts for the full frequency and momentum de- 
pendence of the dynamic screening in the high-density, 
small gradient limit of QCD. The implementation of such 
a self-consistent dynamical screeing mechanism, as for- 
mulated recently by Arnold et al. in a PCM remains 
a major challenge for future work. It should also be noted 
that, strictly speaking, newly produced partons start pro- 
viding screening only after Ar ~ 1 /pT and therefore the 
concept of a screening mass may not be well justified for 
times r < 0. 

In order to determine a lower boundary for the allowed 
range of p™™, we calculate at r = 1/p™™ and plot it 
as a function of p™™ in the bottom frame of figure ||. 
The gray shaded area symbolizes the region in which 
pTOin < and where the procedure to simply cut off 
the interaction is no longer valid. It therefore serves as 
a boundary of the allowed range of p™™ for our calcu- 
lation. We find that the smallest allowed value for p™*" 
is « 0.8 GeV. The inclusion of higher order radiative 
corrections and parton fusion processes into our calcula- 
tion will probably alter this value - we shall study these 
effects in a forthcoming publication [||. 

After having determined the meaningful range of p^'", 
our model is applied to the calculation of the time- 
evolution of the energy-density e generated by hard and 
semi-hard parton-parton interactions at RHIC: the up- 
per frame of figure ^ shows this time-evolution of e for 
different values of p™™, to give an estimate on the sen- 
sitivity of e with respect to p™™. Here, e is calculated 
via 
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FIG. 2: Top: time-evolution of the energy-density for differ- 
ent values of Bottom: rapidity distribution of produced 
partons generated by hard and semi-hard parton-parton in- 
teractions at RHIC. 
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and choosing Rt — 2 fm. The maximum energy density 
obtained in Au-f Au collisions in the PCM approach is 
found to be on the order of 100 GeV/fm'^. A detailed 
analysis shows that for times r < 1.0 fm/c e(r) scales 
with 1/t whereas for later times the scaling function is 
1/t^/'^, most likely signaling a transition from longitu- 



10-' 



10' 

CM 

(5 10^ 
a. 

? 10° 

a. 



10"' 



10"' 



VNI/BMS 



■ (Pt""V=0.5 Gev' 

• (PT""V=1-0Gev' 

♦ (Pt""V=2.0 Gev' 



Qo = (Pt""") 
IVcmI < 0.5 



1 2 3 4 5 6 7 

Pt (GeV) 

FIG. 3: Transverse momentum spectrum of produced partons 
in central Au-|-Au collisions at RHIC for different values of 
^mm NQ^g that the contributions to the spectrum below 
stem from parton rescattering. 



dinal streaming to a three dimensional expansion. One 
might be tempted to attribute this scaling to the onset 
of hydrodynamic expansion. However, the small collision 
rates at later times (see below) preclude this interpreta- 
tion. 

The bottom frame of figure || shows the rapidity dis- 
tribution dN / dy of produced partons for different values 
of 

pmm ^ While the distributions cannot be directly com- 
pared to experimental data, they provide a lower limit 
on the entropy produced via hard and semi-hard inter- 
actions. Since at hadronization entropy can only be pro- 
duced or remain constant, the rapidity density of the 
produced partons provides us with a lower bound for the 
rapidity density of hadrons compatible with the param- 
eters of our calculation. For p™*" = 0.7 GeV the ra- 
pidity density of produced partons at mid-rapidity is ap- 
proximately 650. Assuming parton-hadron duality and 
entropy conservation during hadronization this number 
would translate to approximately 430 charged hadrons. 
In order to be compatible with the data p3[ | , it is impor- 
tant for this number to remain well below the measured 
value, to allow for additional contributions due to intial- 
and final state radiation as well as soft particle produc- 
tion. 

While the absolute values for the energy and rapid- 
ity density show a significant dependence on p™^, the 
overall shape of the transverse momentum distribution 
of produced partons changes only in a subtle manner, as 
can be seen in figure |. The main impact the choice of 
pip™ has, is in the low-px starting point of the respec- 
tive distribution - contributions below p™™ stemming 
from parton rescattering. Naively one would think that 
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the transverse momentum distribution above the chosen 
cut-off should not be affected by the choice of the cut- 
off - however, choosing a low cut-off value increases the 
amount of parton rescattering which visibly modifies the 
shape of the pT spectrum. 
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FIG. 4: Top: collision rates for gluon-gluon scattering as a 
function of time for different values of . Bottom: parton 
collision number distribution for different p™'": the probabil- 
ity for a parton to suffer multiple collisions increases dramat- 
ically with decreasing p™'" . 



In order to quantify the amount of parton scatter- 
ing and rescattering we evaluate the collision rates as 
a function of time and pip™. The upper frame of fig- 
ure ^ shows the collision rates for gluon-gluon scatter- 
ing (glue-glue collisions account for roughly half the to- 
tal number of hard scatterings). The time-evolution of 
the collision rate is clearly peaked at r = fm/c when 
the maximum overlap between the two gold nuclei oc- 
curs. For —0.5 < r < 0.5 fm/c the collision rates show 
a rough symmetry around r = fm/c for large values of 
p™'", indicating the dominance of primary-primary par- 
ton scattering contributions to the collision rate during 
early times if the cross section is small. This symme- 
try gets distorted for decreasing when secondary 
parton rescattering increases. For later times the steep 
decrease in the collision rates flattens out due to parton 
rescattering, exhibiting a far larger sensitivity on p™*" 
(and thus the interaction cross section) than at early 
times. The probability for a parton scattering multiple 
times can be directly evaluated by plotting the parton 
collision number distribution, which can be seen in the 
lower frame of figure |. While for all cases it is most 
likely for a parton to only undergo one hard or semi-hard 
collision, the probability for a parton to suffer multiple 
collisions increases dramatically with decreasing pjp™ . 

In a future analysis, we shall attempt to isolate those 
regions in phase-space in which rescattered partons dom- 
inate and identify experimentally accessible observables 
which are most sensitive to multiple binary perturbative 
rescattering of partons. 

In summary, we have calculated the time-evolution of 
the Debye-screening mass fin for Au-|-Au collisions at 
RHIC in the framework of a lowest-order Parton Cas- 
cade Model. We have used the screening mass to de- 
termine a lower boundary for the allowed range of p™*" 
values. For that range we determined the energy density 
reached through hard and semi-hard processes and found 
that its maximum remains below approx. 70 GeV/fm'^. 
We established a lower bound for the rapidity density of 
charged hadrons which is compatible with the available 
data. The possibility of perturbative rescattering among 
partons has been analyzed and is found to significantly 
increase with decreasing p™"- . 
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